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Abstract 

We follow the consequences of internal equilibrium in non-equilibrium systems that has been 
introduced recently [Phys. Rev. E 81, 051130 (2010)] to obtain the generalization of Maxwell's 
relation and the Clausius-Clapeyron relation that are normally given for equilibrium systems. The 
use of Jacobians allow for a more compact way to address the generalized Maxwell relations; the 
latter are available for any number of internal variables. The Clausius-Clapeyron relation in the 
subspace of observables show not only the non-equilibrium modification but also the modification 
due to internal variables that play a dominant role in glasses. Real systems do not directly turn into 
glasses (GL) that are frozen structures from the supercooled liquid state L; there is an intermediate 
state (gL) where the internal variables are not frozen. Thus, there is no single glass transition. A 
system possess several kinds of glass transitions, some conventional (L— T-gL; gL— t-GL) in which the 
state change continuously and the transition mimics a continuous or second order transition, and 
some apparent (L— T-gL; L^GL) in which the free energies are discontinuous so that the transition 
appears as a zeroth order transition, as discussed in the text. We evaluate the Prigogine-Defay 
ratio n in the subspace of the observables at these transitions. We find that it is normally different 
from 1, except at the conventional transition L^gL, where 11 = 1 regardless of the number of 
internal variables. 



* Electronic address: |pdg@uakron.edu 



1 



I. INTRODUCTION 



A. Previous Results 

In a series of papers, we have begun to develop non-equilibrium thermodynamics starting 
from the second law and ensuring the additivity of entropy as a state function [lifl. The 
central idea in this approach is that of internal equilibrium within a macroscopic system E 
surrounded by an extremely large medium S; the two form an isolated system Sq as shown 
in Fig. [H While the entropy S{t) and the general non-equilibrium thermodynamic potential 



see [2] for more details, such as the non-equilibrium Gibbs free energy G{t) of the 
system exist even when the system is not in internal equilibrium, the Gibbs fundamental 
relation exists only when the system is in internal equilibrium: 

dS{t) = y(t)-dX(t)+a(t)-dI(t), (1) 

where X(t) and l{t) represent the set of observables and the set of internal variables, re- 
spectively, to be collectively denote by Z(t). The entropy S{Z{t),t) away from equilibrium, 
no matter how far from equilibrium, is normally a function of Z(t) and t. However, when 
the system is in internal equilibrium, where Eq. remains valid, S{t) has no explicit 
t-depenedence; the temporal evolution of the entropy in this case comes from the time- 
dependence in Z(t), with X(i(:) and still independent of each other. The coefficient 
y(t) and a(t) represent the derivatives of the entropy and are normally called the internal 
field and the internal affinity, respectively. The energy E, volume V and the number of 
particles play a very special role among the observables, and the corresponding internal 
fields are given by 



j_ ^ f ds{t) \ p{t) ^ fds{ty\ ^ _ ( ds{t) \ 

T{t) \dE{t))^,^,^' Tit) \dV{t)Jz'it)' T{t) \dN{t)J 



(2) 



where Z'(t) denotes all other elements of Z(t) except the one used in the derivative. Thus, 
internal temperature, pressure, etc. have a meaning only when the system comes into 
internal equilibrium. In general, the internal field y{t) and affinity a(t) are given by 

" ^ " (H) ' "^'^ ^W)" im) ■ 

The fields of the medium To,Po,/^o, etc., which we collectively denote by Yq, are different 
from the internal fields of the system unless the latter comes to equilibrium with the medium. 
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FIG. 1: Schematic representation of a system T, and the medium S surrounding it to form an 
isolated system Sq. The medium is described by its fields To,Po, etc. while the system, if in 
internal equilibrium (see text) is characterized by T{t), P{t), etc. 



The same is also true of the affinity, except that the affinity vector Aq = for the medium; 
see II. 

From now on, we will only consider the case when the system is in internal equilibrium. 
The heat transfer is given by 

dQ = T{t)dS{t) = TodeS{t), (4) 

where deS{t) is the entropy exchange with the medium. The irreversible entropy generation 
diS{t) within the system is given by 

diS{t) = dS{t) - d^S{t) > 0. 

Thus, as long as the system is not in equilibrium, T(t) 7^ Tq; accordingly, diS{t) > in 
accordance with the second law. There is irreversible entropy production even when the 
system is in internal equilibrium; the latter only allows us to introduce the internal fields 
and affinities via Eqs. ([2]) and ([3]). 

In the absence of any internal variables, the Gibbs fundamental relation is given by 

T{t)dS{t) = dE{t) + P{t)dV{t) - fi{t)dN{t) (5) 

for the special case when X()f:) only contains E{t),V{t) and N{t). For a fixed number of 
particles, the last term would be absent. As said above, the temperature, pressure, etc. 
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of the medium and the system are usually different when the system is out of equilibrium 
with the medium. Only in equilibrium do they become equal, in which case, the Gibbs 
fundamental relation in Eq. reduces to the standard form 

in which none of the quantities has any time-dependence; the extensive quantities represent 
the equilibrium values and are denoted by the additional suffix. One normally considers a 
system with fixed number of particles, in which case, the last term is absent in Eq. ([6]). In 
the following, we will not explicitly show the additional suffix unless clarity is needed. The 
following Maxwell relations that follow from the Gibbs fundamental relation, see Eq. (E]), 
are well-known and can be found in any good text-book such as on thermodynamics: 

'dTo\ fdPo\ fdTo\ fdV 



9V J g \dSjy^^ \dPoJ \dS / p^ j^ 



\9ToJv,N \9VJto,n \9PoJto,n \9ToJ p,m 



To,N V'^-' U/ To,iV V'^-^U/Po^Tv 

In equilibrium, there is no explicit t-dependence in Z; moreover, the internal variable I is 
no longer independent of X. The equilibrium field and affinity of the system become equal 
to those of the medium (Yq and Aq = 0); see [2]. Thus, the Gibbs fundamental relation 
reduces to 

dS = yo-rfX, (8) 

compare with Eq. (Q. The equilibrium value of the internal variable can be expressed as a 
function of the equilibrium value of X: 

I Icq(-^eq)' 

We now observe the similarity between the Gibbs fundamental relations in Eqs. ([6]) and 
([8]). This strongly suggests that there may also exist analogs of the Maxwell relations or other 
important relations that are based on Eq. ([8]) for a system that, although not in equilibrium 
with the medium, is in internal equilibrium. In this sequel to the earlier papers Im], 
which we denote by I, II, and III, respectively, we develop the consequence of this internal 
equilibrium thermodynamics for important relations such as Maxwell relations, Clausius- 
Clayperon equation, etc. These extensions will play an important role in non-equilibrium 
systems that are nonetheless in internal equilibrium. 



The time-variation of the internal temperature T of a non-equihbrium system such as a 
glass is due to the time dependence of the observable X(t) such as E{t),V{t), etc. and of 
the internal variable I{t). For example, at fixed Tq, the internal temperature will continue 
to change during structural relaxation. The internal temperature will also change if the 
temperature of the medium changes. Thus 

The rate of change of the internal temperature can be expressed in terms of the rate of 
change r = dTo/dt : 

dT 

Similarly, 

dT 

dTo 

The same analysis can be carried out for other internal fields. 
B. Present Goal 

Our aim in this work is to follow the consequences of internal equilibrium in a non- 
equilibrium system to find the generalization of Maxwell's relations, the Clausius-Clapeyron 
relation, and the relations between response functions to non-equilibrium states. We will be 
also be interested in glasses in this work; they are traditionally treated as non-equilibrium 
states. Therefore, we begin with a discussion of what is customarily called a glass and the 
associated glass transition in Sect. HIl A careful discussion shows that the term does not 
refer to one single transition; rather, it can refer to different kinds of transitions, some of 
which appear similar to the conventional transitions in equilibrium, but the other refer to 
apparent transitions where the Gibbs free energy cannot be continuous. There are some 
well-known approximate approaches to glasses. We will briefiy discuss them. We then 
turn to our main goal to extend the Maxwell's relations, where Jacobians are found to be 
quite useful. Therefore, we introduce Jacobians and their various important properties in 
Sect. IIII A[ This is technical section, but we provide most of the required details so that the 
clarity of presentation is not compromised. An important part of this section is to show that 
the Jacobians can be manipulated in a straight forward manner even in a subspace of the 
variables. This is important as the observations require manipulating the observables and 
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dx) ' H ^ [~dT 



dT\ dX fdT 
dx) "dTo^ Vdl 



dl 

di' 

dl 

dTo' 



(9) 
(10) 



not the internal variables. Thus, the experimental space refers to a subspace (Sect. IIIIBp of 
the space where non-equilibrium thermodynamics is developed. Thermodynamic potentials 
for non-equilibrium states are formulated in Sect. IIVI We develop the generalization of 
the Maxwell's relations in Sect. |Vl We discuss generalization of the Clausius-Clapeyron 
relation in Sect. \V1\ where we also discuss the conditions for phase transitions in non- 
equilibrium states. The response functions such as the heat capacities, compressibilities 
and the expansion coefficients and various relations among them are developed for non- 
equihbrium states in Sect. Sect. I VI II The Prigogine-Defay ratio for glasses are evaluated 
at various possible glass transitions in Sect. IVIIII We compare our approach with some of 
the existing approaches in determining the ratio in this section. The last section contains a 
brief summary of our results. 



II. GLASS TRANSITIONS AND APPARENT GLASS TRANSITIONS 

An example of non-equilibrium systems under investigation here is a glass [sl ; see Figs. [2] 
and [31 A supercooled liquid L is a stationary (time- independent) metastable state j^, which 
for our purpose, represents an equilibrium state (by not allowing the crystalline state into 
consideration), and is shown by the curve ABF under isobaric condition at a fixed pressure 
Pq of the medium. We will refer to the equilibrium liquid always as L in the following. In 
contrast, a non-equilibrium liquid state will be designated gL here, and represents a time- 
dependent metastable state [5|. The choice of gL is to remind us that it is a precursor to 
the eventual glass GL at a lower temperature. The equilibrium liquid L is obtained by 
cooling the liquid L and waiting long enough at each for it to come to equilibrium with 
the medium. However, if it is obtained at a fixed cooling rate r, then at some temperature 
Tog(Po)j L cannot come to equilibrium and turns into gL; the resulting curve BD leaves ABF 
tangentially at B, and gradually turns into an isobaric glass GL represented by the segment 
DE at D, when the viscosity becomes so large (~ 10^^ poises) that it appears as a solid. At 
B, the transition is from an equilibrium liquid L to a non-equilibrium liquid form gL, and 
will be called the L-gL transition. In the literature, it is commonly known as a transition 
from an ergodic state (L) to a non-ergodic state (gL). In our opinion, this is a misnomer, as 
the concept of ergodicity refers to the long-time , indeed the infinite-time, behavior. In this 
limit, there will be no gL, only L. Therefore, we will refer to this transition at T = Tog(Po) 
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as a L-gL transition or a precursory glass transition. The true glass transition at D is not a 
transition from L to GL, but a transition from gL to GL. We will refer to the glass transition 
at the lower temperature Tog(-Po) at D as the actual glass transition, or simply the glass 
transition. The transition region BD represents a time-dependent metastable supercooled 
liquid (to be distinguished from the stationary metastable supercooled liquid L denoted by 
ABF), which turns into a glass at D. The expansion coefficient in the glass is almost identical 
to that of the corresponding crystal below D. The glass continuously emerges out of gL at 
D, whose location is also determined by the rate r of cooling. The relaxation time r of 
the system (the supercooled liquid) becomes equal (really comparable) to the observation 
time Tobs at B. As seen in Fig. m the volume remains continuous at B and D at the two 
glass transitions. The same is also true of the entropy. Indeed, the state of the system 
changes continuously at B and D, which is highly reproducible for a given cooling rate r 
or the observation time Tobs- Thus, the points B and D can be taken as a well-defined 
and unique glass transition temperatures Tog(Po) and Tog(-Po) associated with the point B 
and D, respectively, in both figures. Both transitions represent a non-equilibrium version 
of a continuous transition (See Sect. |Vl]for elaboration on this point), where not only the 
Gibbs free energy, see Fig. [3l but also its derivatives are continuous. The non-equilibrium 
nature of the transition appears in the dependence of the value of Tog(Po) and Tog(-Po) on 
the rate of cooling. The continuity of the Gibbs free energy at B and D makes them as 
genuine candidates as (glass) transition points, a requirement of a transition in equilibrium 
thermodynamics. Therefore, both these transitions will be collectively called conventional 
transitions in this work. 

Unfortunately, the idea of a glass transition was formulated as a transition between L 
and GL. Thus, neither of the above two glass transitions represent the glass transition in 
the original sense. As the glass is considered a frozen state, it is common to assume that 
over the region DE, the glass has its internal variables denoted by I frozen at its value Iq at 
D, even though its observables denoted by X continue to change. On the other hand, the 
internal variables and the observables continue to change over BD from their values at B to 
their values at D. Consequently, the properties such as the volume of gL, which is shown 
schematically in Fig. [2], gradually change to those of the glass at lower temperatures. Thus, 
the glass transition from AB to DE is not a sharp transition. It can be argued, as we have 
done above, that B and D should be taken as the glass transition points. However, the 
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FIG. 2: Schematic form of isobaric V as a function of Tq for a given cooling rate. The pressure is 
fixed at Pq. The supercooled liquid turns gradually into a glass through the glass transition region. 
The transition temperature Tog(Po) is identified as the temperature at B, where the actual volume 
begins to deviate from the extrapolated supercooled liquid volume BC. On the other hand, the 
apparent glass transition temperature Tq^\Pq) is the temperature where the extrapolated glass 
volume DC meets the extrapolated supercooled liquid volume BC as indicated in the fuigure; this 
temperature lies in the glass transition region. 



practice in the field is to take a point between BD as a transition point obtained by electing 
some well-defined rule of selection; see for example |6] for a good discussion of various ways 
of identifying the glass transition temperature. One such rule commonly used is to consider 
the volume of the system and introduce an apparent glass transition temperature Tq^\Pq) 
by the equilibrium continuation of the volume BCF of AB and by the extrapolation of the 
volume DCG of DE to find their crossing point C. The state of the glass following Tool 
and Narayanaswamy 0, sl is then customarily identified by the point C on DC. However, 
there is no reason to take the state at C to represent any real glass, as the extrapolation 
does not have to satisfy non-equilibrium thermodynamics; the latter is valid only along the 
physical path DB for the given history of preparation such as determined by the fixed rate 
r of cooling during vitrification. The glass at Tq^\Pq) must be described by the point on 
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FIG. 3: Schematic form of the isobaric Gibbs free energy G shown by the continuous curve ABDE 
as a function of the medium temperature Tq at a fixed pressure Pq. The extrapolation of the glassy 
portion (GL) along DCG and the supercooled liquid (L) portion ABCqF do not meet; the glassy 
Gibbs free energy at point the apparent glass transition C, where Tq = Tq^\Pq), is higher than that 
at Co on the continuous curve L at the same temperature Tq^\Pq), showing that the extrapolation 
results in a more unstable state at the apparent glass transition C than the physical state Co on 
the continuous curve. The Gibbs free energies match at the glass transition temperature rog(Po) 
at B. 

DB corresponding to Tq^\Pq) if we wish to employ non-equilibrium thermodynamics. To 
be sure, one can find a slow enough cooling rate than the one used to obtain gL at B so 
that the point B actually coincides with the point C on ABF, as the latter represents L. 
However, the gL that will emerge at C for the slower cooling rate has nothing to do with the 
extrapolated state C on DCG. Because of the continuity of the state, the gL at the slower 
rate at C will have its A = and ^ = ^eq and will have its Gibbs free energy continuous. 
Moreover, the new gL will follow a curve that will be strictly below BDE. These aspects 
make the new gL different from the extrapolated GL at C. Taking the point C on CD to 
represent the glass will be an approximation, which we will avoid in this work, as our interest 
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is to apply thermodynamics in the study of glasses. Therefore, we will use the extrapolation 
to only determine Tq^\Po), but the real glass and the real liquid states are determined by 
the curve BD and BCF, respectively, where our non-equilibrium thermodynamics should be 
applicable. 

The location of this temperature Tq^\Po) depends on the property being extrapolated. 
We can use the entropy of the system to locate the apparent glass transition temperature, 
which would invariably give a different value for the apparent glass transition temperature. 
To call one of these temperatures as a transition temperature is a misnomer for another 
reason. None of these temperatures represent a " non-equilibrium" thermodynamic transition 
for the simple reason that the two branches DCG and BCF do not have a common Gibbs 
free energy at Tq^\Pq) as is clearly seen in Fig. [31 The branch ABCqF represents the 
Gibbs free energy of the equilibrium supercooled liquid, while the segment DE represents 
the Gibbs free energy of the glass, with the segment BD denoting the Gibbs free energy of 
the system during the transition region. The extrapolation DCG in Fig. [2] to determine the 
glass transition temperature Tq^\Po) corresponds to the extrapolated segment DCG in Fig. 
[3l The Gibbs free energy of the glass in this extrapolation is given by the point C, while the 
Gibbs free energy of the supercooled liquid is determined by the point Cq. Evidently, the 
two free energies are very different, with that of the glass higher than that of the supercooled 
liquid, as expected from the non-equilibrium nature of the glassy state. 

The above discussion of the apparent glass transition also applies to comparing the glass 
at D with the corresponding L at Tqq^Pq), which will represent yet another apparent glass 
transition temperature. This apparent glass transition has the same problem regarding the 
Gibbs free energy as the previous one at Tq^\Pq). However, this transition differs from 
the apparent glass transition at T^^\Pq) in that the "glass" at Tq^\Po) is not a frozen 
state, while the glass at D is a "frozen" glass to some extent (as it also undergoes structural 
relaxation in time). It should also be remarked that whether we consider the apparent glass 
transition at Tog(-Po) or Tq^\Pq), the transition is an example of a discontinuity in the Gibbs 
free energy of the two states. This is different from the precursory glass transition and the 
actual glass transition at B and D, respectively, where the Gibbs free energy is continuous. 
Because of the discontinuity in the Gibbs free energies in the apparent glass transitions 
at Tog(-Po) and Tq^\Po), we will refer to these transitions as apparent transitions in this 
work. Indeed, one can think of these transitions as an analog of a zeroth order transition 
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because of the discontinuity in the Gibbs free energy. However, it should be remarked that 
the apparent transitions do not represent any transition in the system; those transitions are 
the two conventional transitions discussed above. The apparent transitions represent our 
desire to compare two distinct states. This is like comparing the supercooled liquid with the 
crystal at the same temperature and pressure. Therefore, a discontinuity in the Gibbs free 
energy is not a violation of the principle of continuity discussed in j^. 



We will consider all of the a 



of the Prigogine-Defay ratio 



30ve glass transitions later when we discuss the evaluation 



13( 1 in Sect. IVIIII In this ratio, a non-equilibrium state is 



compared with the equilibrium supercooled liquid state along ABF. In the classic approach 
adopted by Simon Q, the temperature range (TogG(-Po)) ^og(-Po)) is shrunk to a point, 
either by considering the apparent glass transition at Tq^\Pq), or by comparing the glass 
state at D with the supercooled liquid L at B. The latter amounts to neglecting the segment 
BD from consideration. We will avoid this ad-hoc approach in this work. The only possible 
scenario, where Simon's approach is meaningful is that of the ideal glass transition and 
references thererin], in the limit the cooling rate r — )■ 0. In this limiting case, the crossover 
region BD disappears and the ideal glass IGL emerges directly out of the L at the ideal 
glass transition temperature Tqig- This is a conventional continuous transition between the 
two stationary states IGL and L, both of which remain in equilibrium with the medium at 
To, Pq- There is no need to invoke any internal variable I to describe the ideal glass; the 
observable X is sufficient for the investigation of the ideal glass transition. We will revisit 
this point later in Sect. IVIIII 



III. SOME USEFUL MATHEMATICAL TOOLS 



Jacobian method 



Jacobians 15| will be found extremely use: 
equilibrium thermodynamics ^; see also 
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'ul in this work just as they are found useful in 
-|l8|. The n-th order Jacobian of Ui, U2, ■ ■ - Un 
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with respect to xi,X2, ■ • • Xn is the n x n determinant of the matrix formed by duk/dxf. 



d{Ui,U2, ■■■Un) 
d{xi,X2, ■■■Xn) 



dui/dxi dui/dx2 
du2/dxi du2/dx2 



dui/dxn 

du2/dXn 



dUn/dXi dUn/dx2 ■ ■ dUn/dXn 

It is clear from the properties of the determinant that 
1. The Jacobian vanishes if any two tt's are identical 

d{ui. 112- ■ ■ ■ ti,. u, ■ ■ ■ (/„) 



, X2i ■ ■ ■ Xi^ ^-i+l ■ ■ ■ ^^n) 



0. 



2. If ti j and tij+i interchange their order, the Jacobian changes its sign 

d{ui,U2, ■ ■ -Uj+i^Uj ■■•Un) _ d{ui,U2, ■ ■■Uj^Uj+i ■ ■ ■ Un) 
1 X2i ■ ■ ■ Xj^i • • • Xn) di^Xi, X2i ' ' ' X^^ ■ ■ ■ X^) 

3. If any Ui is equal to Xj, the n-th order Jacobian reduces to a (n — l)-th order Jacobian 
formed by derivatives at fixed Xj. For example, for n = 2, we have 

d{ui. .V2) ( du\ 



d{xx,X2) \dxx 

When we consider compound transformations (xi, 2:2, • • -x^) — >■ (i^i, ^2, • " •'^n) 
{u\^V2i^ ■ ■ Vn)i the resulting Jacobian is the product of the two Jacobians: 

d{yx,V2-, ■■■Vn) _ d{Ui,U2, ■■■Un) _ d{Vi,V2, ■■■Vn) 
d{ui,U2r • - Un) d{xi,X2,^ ■ ■Xn) d{xi, X2, ■ ' ' ^n)' 

The definition of a Jacobian can lead to some interesting permutation rules as the 
following examples illustrate. Consider a second order Jacobian d{ui,U2)/d{xi,X2) — 
(dui/dxi) {du2/dx2) — {dui/dx2) {du2/dxi), which can be rearranged as 

d{ui,U2) d{xi,X2) , d{u2,Xi) d{ui,X2) , d{xi,Ui) d{u2,X2) 



+ 



+ 



9(xi,X2) 9(xi,a;2) 2:2) 3:2) 9(a;i,a:;2) 9(a:;i,a;2) 

This can be symbolically written as 



= 0. 



d{ui,U2){xi, X2) + d{u2, xi)d{ui,X2) + d{xi,ui)d{u2, X2) = 



(11) 
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by suppressing the common denominator in each term. The result expresses the cychc 
permutation of ui, U2, xi in the three terms with the remaining variable X2 in the same place 
in all terms. As a second example, consider some quantity m as a function of three variables 
X, y, and z and consider the following relation between the partial derivatives: 

(^) , ^ (^) ^ (^) (^) , ■ ^^^^ 

In terms of Jacobians, it can be written as 

d{u,y) ^ d{u,y,z) ^ d{u,x,y) d{z,y) 
d{x, y) d{x, y, z) d{z, x, y) d{x, y) ' 

which simplifies to 

d{x, y, z)d{u, y) = d{y, z, u)d{x, y) + d{z, u, x)d{y, y) + d{u, x, y)d{z, y), (14) 

where we have added a vanishing second term on the right because d{y, y) = 0. This relation 
is easily constructed by considering the cyclic permutation of 

x,y,z,u 

by taking three consecutive terms at a time for the 3- Jacobians, with the remaining variable 
yielding the 2- Jacobians in which the second entry is the variable y, the variable that is 
held fixed in all derivatives in Eq. f[T^ . The ordering x,y,z in x,y,z,u is determined by 
the denominator 3-Jacobian in the first term on the right in Eq. f|T3l) . By writing all the 
3- Jacobians in the non-vanishing terms in Eq. (IT^ so that y is the second entry, and then 
suppressing the second entry, we obtain the following relation 

d{x, z)d{u, y) + (9(z, u)d{x, y) + d{u, x)d{z, y) = 0, 

which is identical to the relation in Eq. ( ITT]) if we identify ui with x, U2 with z, xi with u 
and X2 with y. 

We will use the Jacobians and their properties to first re-express the Maxwell relations 
as follows 

9(To,5,iV) d{P^,V,N) d{To,S,N) d{Po,V,N) 



d{V, S, N) d{V, S, N) ' d{Po, S, N) d{Po, S, N) ' 
d{Po,V,N) _ diTo,S,N) d{To,S,N) _ diPo,V,N) 
d{To,V,N) diTo,V,Ny d{Po,To,N) 9(Po,To,iV)' 



(15) 
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We now see a very important consequence of the use of the Jacobians. All four Maxwell 
relations use the same numerators d{To, S, N) and d{Po, V, N). They use different denomi- 
nators. Thus, they can all be combined into one compact relation that can be simply written 

as 

d{To,S,N) = d{Po,V,N). (16) 

Here, the relation only has a meaning if each side is divided by one of the possible denomi- 
nators d{V, S, N), d{Po, S, N), d{To, V, N) and d{Po, To, N) on both sides. 

B. Considerations in a Subspace 

It is very common to consider a function F{x, y, z) in a subspace consisting of x, y. This 
requires manipulating a 3- Jacobians to construct a 2- Jacobians. of its argument. Thus, we 
may consider the 2-Jacobian 

d{x, y) ' 

even though F also depends on z. We can manipulate such Jacobians in the normal way. 
For example, we can express it as 

(dF\ _ d{F,y) _ d{F,y)d{x,K) _ fdK\ d{F,y) 

\dxj^ dix,y) diK,x)dix,y) \dy ) ^d{K,xy ^ > 

where K[x,y,z) is another function. The derivation is tedious and has been supplied in 
the Appendix. The situation can be generalized to many variables zi,Z2, - ■ ■ without much 
complications. We will not do this here. 

C. Some Transformation Rules 

Let us consider a derivative of some quantity R either with respect to T or P in case A 
below or at fixed T or P in case B below, which we wish to express as a derivative involving 
To, Pq that are manipulated by the observer. 

A. The derivative is at fixed U, where U has any two different elements from E,V, S,^, Pq 
and To- 
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We write the derivative as 



(18) 

Similarly, we have 



9PJu,,u, ~ d(P,U,,U2) d{Po,U,,U2) d{P,U,,U2) \dPj^^ J \dPj 



Ul,U2 

(19) 



B. Let us consider a derivative with respect to Tq at fixed U2 — P or T {U20 — Pq or Tq, 
as the case may be), but Ui is any element from E, V, 5", ^, Pq and Tq: 

dR\ _ d{R,U,,U2) d{To,U,,U2o) ^ d{R,U2,Ui) , ( dU2\ , . 

^oJu,,u, ~ 9(To, C/i, U20) d(To, U,, U2) d(To, U20, U,)' \dU2oJu,,u, ' 

Similarly, 

97? \ _ diR,UuU2) d{Po,U^M2o) _ diR,U2,Ui) J dU2\ . . 

dPoJu„u,~ 9iPo,U,,U2o) d{Po,U,,U2) d{Po,U2o,U,)^ \dU20 J u„u,' 

Let us now consider a derivative with respect to T at fixed P or with respect to P at 

fixed T; the derivative is at fixed C/i, where Ui is any element from E,V, S,^, Pq and 

T 

'dR\ _ d{R,U,,P) _ diR,U,,P) diT,U,,P) 
dTj^^p d{T,Ui,P) d{To,U,,Po)^ d{To,U,,Poy ^ ^ 



dR\ _ d{R,UuT) d{{R,Ui,T) d{P,Ui,T) ^^3) 



Similarly, 

{dPj^^,.- d{P,U^,T) d{To,U^,Po)' d{To,U,,Po) 
IV. THERMODYNAMIC POTENTIALS AND DIFFERENTIALS 



A. Equilibrium 

The forms of most useful thermodynamic potentials such as the enthalpy H, the 
Helmholtz free energy F, and the Gibbs free energy G of a system E in equilibrium are 
well known and are given in terms of the energy E{S,V,N) as 

H = E + PoV, F = E- ToS, G = E-ToS + PqV, (24) 

where Tq, Pq are the temperature and pressure of the system; they are also the temperature 
and/or pressure of the medium, depending on the medium E. Here, we are considering 
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a system with fixed number of particles. For tlie entfialpy, tlie medium S(Po) containing 
tlie system exerts a fixed pressure Pq. For the Helmholtz free energy, the medium S(To) 
containing the system creates a fixed temperature Tq. For the Gibbs free energy, the medium 
S(To,Po) containing the system exerts a fixed pressure Pq and creates a fixed temperature 
Tq. The potentials are Legendre transforms in that the potentials are functions of the fields 
(Tq, Pq) rather than the observables {E, V) as the case may be. These potentials have the 
desired property that they attain their minimum when the system is in equilibrium, as 
discussed in I. 

B. Internal Equilibrium 

When the system is in internal equilibrium, we find from the Gibbs fundamental relation 
for fixed A^, which is obtained from setting dN = in Eq. (|5]): 

dE = TdS - PdV - Adi, (25) 

where we have also introduced a single internal variable ^ to allow us to discuss non- 
equilibrium systems that are not in equilibrium with their medium but are in internal 
equilibrium. The consideration of many internal variables is to simply replace 

Adi A-dC,, 

and will not cause any extra complication. Thus, we will mostly consider a single internal 
variable, but the extension to many internal variables is trivial. 

We are no longer going to exhibit the time-dependence in these variables for the sake 
of notational simplicity of. Let us return to Eq. ( I25l) . It should be compared with Eq. 
dH]) which contains Tq, Pq. We rewrite Eq. (1251) to show the non-equilibrium contribution 
explicitly: 

dE = TodS - PodV + (T - To)dS - (P - Po)dV - Ad^. (26) 

The last two terms are due to the non-equilibrium nature of the system in internal equilib- 
rium. It is now easy to see that 

dH = TodS + VdPo + (T - To)dS - (P - Po)dV - Ad^, 

dF = -SdTo - PodV + (T - To)dS - (P - Po)dV - Ad^, (27) 

dG = -SdTo + VdPo + (T - To)dS - (P - Po)dV - Ad^. 
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These potentials correspond to ^ as an independent variable of the potential. One can make 
a transformation of these potentials to potentials in which the conjugate field Aq of the 
medium is the independent variable by adding Aq^. The resulting potentials will be denoted 
by a superscript A on the potential: 

E^ = E + Ao^, H^ = H + Aoi, = F + Aq^, = G + Aq^. 

However, as discussed in II, Aq = 0. Thus, there is no difference in the values of the two 
potentials and the transformation is of no use. In equilibrium, the internal fields T, P attain 
their equilibrium values Tq, Pq of the medium, and the affinity A vanishes identically because 
of Aq = 0. 

V. MAXWELL RELATIONS FOR SYSTEMS IN INTERNAL EQUILIBRIUM 

From now on, we will always consider the case of a constant A^. Therefore, we will no 
longer exhibit it anymore. The Maxwell relation in Eq. f lT6|) will then be denoted simply 
as d{TQ,S) = d{PQ,V). The field parameters that appear in the Maxwell relation are the 
parameters Tq, Pq of the medium, which because of the existence of equilibrium also represent 
the field parameters of the system. The Maxwell relation is a relation between the pairs 
To, S and Pq, V, each pair formed by the extensive variable and its conjugate field. We will 
call these pairs conjugate pairs in this work. For a system described by only two conjugate 
pairs, there is only one possible Maxwell relation. For a system described by three conjugate 
pairs, there will be three different Maxwell relations between them. For a system described 
by k conjugate pairs, there will be k{k — l)/2 different Maxwell relations. 

As the system in internal equilibrium is very similar in many respects with an equilibrium 
system as discussed in I and II, there may be analogs of the Maxwell relations for systems in 
internal equilibrium. The question then arises as to the field parameters that must appear in 
the Maxwell relations when the system is not in equilibrium, but only in internal equilibrium. 
We now turn to answer this question. Because of the absence of equilibrium, we must now 
also include the internal variable ^ in the discussion. Thus, we expect three different Maxwell 
relations between 
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A. Maxwell relation d{T, S, C) = d{P, V, C) at fixed C 



We start with Eq. (125|) and observe that 



Using the first two derivative at fixed ^, we find that 



\dVds)~\dv)s.'\dSdv)~ \dS ) . 



As we are allowed to interchange the order of derivatives in the above cross derivative, we 
have 

which can be written using Jacobians as 

d{T,S,0 _ diP,V,0 

d{s, V, d{s, V, ■ 

This suggests the existence of the Maxwell relation d{T,S,C,) = d{P,V,^) between the 
conjugate pairs T, S and P, V at fixed C,- To check its validity for other potentials with ^ as 
an independent variable, we consider the differential dG in Eq. (1271) and note that 



We use these derivatives to evaluate the cross derivative (d'^G/dPodTo)^ to conclude that 



( d^v \ (dv\ _ fdP_\ fdv_ 

( d^V \ _ fdS_\ _ fdP\ fdV\ 
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This is simplified to yield 



In terms of the Jacobians, this can be written as 



thus justifying the Maxwell relation 

d{T,S,0=d{P,V,0 (29) 

at fixed ^. This relation must be satisfied at every point on the curve ABDE that describes 
the vitrification process. This Maxwell relation turns into the identity 



d{T,S,0 _ diP,V,0 

for the enthalpy and 

d{T,S,0 _ d{P,V,0 
diTo,V,0 d{To,V,0 
for the Helmholtz free energy, and are easily verified. 



(30) 



B. Maxwell relation d{T, S, V) = d{A, ^, V) at fixed V 

We again start with Eqs. and fl25]) , and evaluate the cross derivative {d'^E /dSd^)y 
to obtain 

/ d^E \ _ fdT\ / d'^E \ _ _ fdA\ 

We thus have 

dT\ _ _ fdA\ 
sy J v,i 



which can be written using Jacobians as 

d{T,S,V) d{A,i,V) 



19 



This suggests the existence of the Maxwell relation d(T, S, V) = d{A, ^, V) between the 
conjugate pairs T, 5* and A, at fixed V. To check its validity for other potentials with V 
as an independent variable, we consider the differential dF in Eq. (1271) and note that 



dF 



-S+{T-To) 



0/ v,^ 



dS 



fdF_ 



0/ V,i 



We now evaluate the cross derivative [S^F /d^dTQ)y and obtain the equality 

dS 



+ 



Toy 
OA 



+ 



0/ v,^ 

which leads to the relation 



+ {T-To: 



To,V 



Wot J. 



'4) HT-T^ 



To,V 



diA,^,V) diT,S,V) 



\d^dTo 



0/ V 



This confirms that the Maxwell relation between the conjugate pairs T, S and A,C, at fixed 
V is the following: 

d{T,S,V)=diA,^,V). (31) 



C. Maxwell Relation d{P, V, S) = d{A, ^, S) at fixed S 

We again start with Eqs. ( 125|) and (!28l) . and evaluate the cross derivative {d'^E/dVd^)g 
to obtain 



We thus have 



dP\ _ rdA 

which can be written using Jacobians as 

d{P,V,S) _ diA,^,S) 
d{^,V,S) d{^,v,s)- 

This suggests the existence of the Maxwell relation d{P, V, S) = —d{A, ^, S) between the 
conjugate pairs P, V and A, ^ at fixed S. To check its validity for other potentials with 5" 
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as an independent variable, we consider the differential dH in Eq. ( 127|) and note that 



V^/ 5.P0 



We now evaluate the cross derivative {d^H/dC,dPo)g and obtain the equality 

- 1 



OA 

which leads to the relation 



dP 



dPo 



/dV_ 

V^/ Po,S 



d{A,^,S) d{P,V,S) 



d{Po,^,s) d{Po,^,sy 

This confirms that the Maxwell relation between the conjugate pairs P, V and A, ^ at fixed 
S is the following: 

diP,V,S) = -diA,^,S). (32) 

One can easily check that this Maxwell relation also works with other thermodynamic 
potentials like F and G. We will satisfy ourselves by giving the demonstration for F only. 
The natural variables for F are Tq, ^ and V; however, instead of using V as the independent 
variable, we will use S as the independent variable so that it can be held fixed. For constant 
S, the differential dF from Eq. ( 127|) reduces to 



dF\ 



-SdTo - PdV - Adi, 



so that 



dF 



-S-P 



0/ 5,C 



dV 



Now evaluating the cross derivative {d'^F/d^dTQ)g, we find that 



dA\ fdP\ fdV_\ d^V \ 



dr. 



dV\ /dP\ 



To,S 



\d^dPo), 
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This now immediately leads to 

d{A,^,S) ^ d{P, V, S) 

d{To,^,S) d{To,^,Sy ^ ' 

and confirms our claim that the Maxwell relation is given by Eq. (!32|) . 

D. General Maxwell Relations with system variables only 

We wish to emphasize that the Maxwell relation in Eq. f p3|l requires keeping S fixed so 
that we must divide Eq. ( 132|) by 9(To, ^, S) on both sides. We must not use the independent 
variables Tq, ^ and V of F for the division and keep To fixed. This will not give be a Maxwell 
relation. We demonstrate this explicitly by evaluating {d'^F/dC,dV)rp^^ two different ways and 
equating the results. A simple calculation yields 

V,To V^'' / To,5 V / V,To V'^'' ^"^y To 



fdS\ fdT\ ,^ f d^S \ 
faA\ fdS\ (dT\ ^\f<>''S\ 



In terms of Jacobians, the above equation can be rewritten as 

d{A,^,To) ^ diP,V,n) diT,S,n) 

d{V,^,To) d{V,^,To) d{V,^,To)- ^ ^ 

This relation from the cross derivative requires keeping Tq fixed. However, Tq characterizes 
the medium and only indirectly characterizes the system in internal equilibrium. In a similar 
way, using the cross derivatives of the Gibbs free energy at fixed Tq, and at fixed Pq, we find 
the following relations: 

d{A,^,Po) d{T,S,Po) d{P,V,Po) d{A,^,To) d{P,V,To) d{T,S,To) 
d{To,Po,0 5(To,Po,0 5(To,Po,0' '9(To,Po,0 5(To,Po,0 '9(To,Po,0' ' 

We now wish to observe that the Maxwell relations appear only when we keep the quantities 
of the system T, P, S, V, A, or ^ fixed. We have already seen the Maxwell relations with fixed 
S, V, and ^. We will now consider keeping T fixed to demonstrate our point. For fixed T, 
we obtain the following Maxwell relation 

diP,V,T) 



5(To,e,T) 9(To,e,T)' 



(36) 
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as can easily be checked by evaluating the cross derivative {d'^F/d^dTo)j, at fixed T. The 
calculation is identical to that carried out in obtaining Eq. (|33|) . One can easily check that 
keeping P or A also gives us new Maxwell relations 

d{A^,P) ^ diT,S,P) diT,S,A) ^ d{P,V,A) 

9(To, p) 9(To, e, p) ' To, A) din, e, a) • 

VI. CLAUSIUS-CLAPEYRON RELATION 

As a system in internal equilibrium is not very different from that in equilibrium, except 
that its Gibbs free energy G{t) continuously decreases until it reaches equilibrium with the 
medium, it is possible for the system to exist in two distinct phases that have the same 
Gibbs free energy at some instant. Such a non-equilibrium phase transition situation will 
arise, for example, when an isotropic supercooled liquid can turn into a liquid crystal phase. 
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, and references 



This is not a novel idea as there are several attempts in the literature 
therin] where such non-equilibrium phase transitions have been investigated. Therefore, let 
us now consider the possibility of the system being in two different phases at some time. As 
experiments are carried out by controlling observables only and not the internal variables, it 
is important to consider thermodynamic quantities as a function of X only, and not of X, I in 
all cases. Restricting ourselves to a single internal variable ^, and to E and V, we will treat 
thermodynamic quantities not only as a function of three independent variables, but will 
also have the need to consider them as a function of observables or associated fields To, Pq- 
In particular, the Clausius-Clapeyron relation is obtained in the Tq-Pq plane, a subspace; 
see Sect. UTTBI 

Let us consider the two phases, which we denote by 1 and 2, in the system. We will 
use subscripts 1 and 2 to refer to the quantities in the two phases. In internal equilibrium, 
the entropy S of the system is a function of the averages X(t),I(t) along with the fixed 
number of particles A^. It is important to include in our consideration as the two phases 
will contain number of particles Ni and N2 that are not constant, except in equilibrium. 
Obviously, 

X(t) = Xi(t) + X2(t), I{t) = hit) + Ut),N = N,{t) + N2{t). 
Then, we can express the entropy of the system as a sum over the two phases: 
S{X{t), I{t), N) = S,{X,{t), hit), N,{t)) + S2{X2{t)Mt), Mt)), 
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which takes its maximum possible value for given X(t), I(t), in internal equilibrium. Thus, 

dS{X{t), I(t), N) = dS^iX,{t), hit), iVi(t)) + dS2{X2it), h{t), N^it)) = 

in internal equilibrium. This can only happen if 

yi(t) = y2(t), ai(t) = a2(t), /ii(t)/Ti(t) = ^2(t)/T2(t); 

see Eqs. (HD, © and (E]). 

For the restricted case under consideration, this results in the equality 

Ti(t) = T2(t), Pi{t) =P2(t), l^i{t) = /i2(t), Aiit) = A2{t) 

for the internal fields and affinity along the coexistence of the two phases. It also follows 

n 

from the continuity of the Gibbs free energy p| that the Gibbs free energies of the two pure 
phases (A^^i = and N2 = N) must be equal at the coexistence. We will only consider the 
two pure phases below, and not a mixture of the two. As the numbers of particles in the 
two pure phases are constant, we will no longer consider them anymore in the discussion. 

We now consider the Tq-Pq plane, relevant for the observation of coexistence. Since the 
Gibbs free energy is continuous along the transition line, 

AG'(To,Po(To)) = 
where Po(To) is the pressure along the transition line. Thus, 



Using dAG = yields 



coex 



dPo 

along the coexistence. Using dG from Eq. (!27|) gives us 



(39) 

To 
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Putting the above two equations in Eq. (1371) . we get the following Clausius-Clapeyron 
equation for coexistence of phases in internal equilibrium 



dTn 



AV + {T- To) A {dS/dPo)T, -{P- Po) A {dV/dPo)^^^ - AA {d^/dPo)^^ 



AS-{T- To) A idS/dTo)p^ + (P - Po) A {dV/dTo)p^ + AA {d^/dTo)p^ ' 

(40) 

We now express {dS/dPo)rp^ in terms of {dV/dTQ)p^ by using the Maxwell relation d{P, V) = 
d(T, S) and by using Eq. ( fTTl) {F S, K ^ V,x ^ Pq, and y — )■ To) as follows: 

d{S,To) ^ d{S,To) d{P,V) d{Po,V) 

9(Po, To) diS, T) diPo, V) diPo, To) ' 

which immediately gives 

ds\ {dP/dPo)y (dv \ ^^^^ 



which can now be used in the Clausius-Clapeyron equation to express it in terms of mea- 
surable quantities assuming that P, T can be measured. In equilibrium, T = Tq, P = Pq and 
A = 0, so that the above equation reduces to the well-known version 

(cq) 



dTn 



AS' 



dPo 

as expected. 



VII. RESPONSE FUNCTIONS IN INTERNAL EQUILIBRIUM 

A. Cp and Cy 

The heat capacities with respect to the internal temperature at fixed P or \^ are 
We again start from the fundamental relation in Eq. (jlHjl and evaluate the derivative 



(42) 



dTjp^ \dT)y^^ \dV)^^\dT)p^ 
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which can be rewritten in two equivalent forms 
C = Cv,g + T 



dS 
dP 



df) 



I 



or 



C 



C 



T 



\dT Jy^^ \dT 



(43) 



(44) 



where we have used the Maxwell relation in Eq. fl29|) after we divide it by d{V,T,C,)- As 
{dS/dP)y^ is not directly measurable, the identity in Eq. (jH]) is more useful from a practical 
point of view. However, we need to transform the various derivatives in it to the derivatives 
with respect to Tq at fixed Pq V by using the transformation rules in Sect. IIII CI as it 
is the pair Tq,Pq that can be manipulated by the observer. However, the identities still 
contains C'p,^ and 6*^,^, which are defined with respect to T, and not with respect to Tq. 
Therefore, we now turn to heat capacities obtained as a derivative with respect to Tq. 



B. Cp and Cy 



From Eq. (jll), we have 



c.--m --rm .Tf#i , (45) 



It would have been more appropriate to express the capacities Cp and Cp^^ as Cpg and 
Cpg ,c, but we will use the simpler notation. This should cause no confusion. Introducing 
the expansion coefficient 

1 fdV\ 1 fdV\ 

. api: = — —— 



we find that 



Cp ^ d{S,Po)/d{T,,Po) ^ d{S,Po) ^ fdS_ 
ap d{V,Po)/d{To,Po) d{V,Po) \dV 



Po 



The same discussion can be applied to Cp^^ and ap^^ with a similar result 

Cp^ ^ d{S,Po,0/d{To,Po,0 ^ d{S,Po,0 ^ fdS\ 
ap,^ d{V,Po,0/d{To,Po,0 d{V,Po,0 \9V J p^^^ 
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Let us now consider the relation between Cp^^ and Cy,^ and between Cp and Cy, for 
which we consider S* as a function of T, V and S,-, which follows from Eq. fl25|) . so that 

= + ai^^^ + ^ ^ 

Therefore 

_ fdS_\ f&T\ fdS_\ fdV\ 
Now, using Eq. (ITSl) . we have 

Similarly, using the Maxwell relation in Eq. ( I30ll we have 



9Ty,., V9Tny,./U^oJ. ^^^^ 



fdS\ ^ d{V,P,0 d{V,To,0 ^ fdP\ fdT\ 
\dVj^^^ d{V,To,0 d{V,T,0 {dToJv,/ WoJv,i' 



We thus finally obtain 



dS\ f dT\ f dS\ f dT\ f dP\ f dV 

+ 



dT^J p^^^ ydToJ ydToJ \dToJ p^^^ \&TqJ y^^ \dTo/ ^ 

After multiplying by T on both sides, we obtain the desired relation between Cp_^ and Cy,^ 
for the non-equilibrium case 

(§) = ^^''^ (^) p./ K§) V. (i^ ) P.. ■ 

This relation generalize the following standard equilibrium relation: 



dToJv \dTo 



obtained by setting 

dT\ f dT 



V \^"U/ Po 



1. 



dT^Jy^^ V^^o/Po,^ 

We can obtain a standard form of the above heat capacity relations as in Eq. (H3|) as 
follows: 

_ d{S,V,i) _ d{S,V,0/d{To,Po,0 
a(To,V^,0 a(To,V^,0/5(To,Po,0' 

We thus finally have 
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which is an extension of Eq. ( H3|) . Although tedious, it is straightforward to show that this 
relation is identical to the above relation. One needs to evaluate {dS/dV)rp^^ as follows: 

9(5, To, 9(V^,T,0 



dS_\ 

dv)^^^^ d{V,T,0 d{V,n,0 



dTo) 



dTo\ 



_ fdS 



where we must now use Eqs. (150|) and (149|) . We finally obtain 



V^/ To,, 



dT 



0/ v,i 



dP 



dTn 



on 



dS 



dTo 



dT0jy^^\dT Jy^^\dV 



dTn 



(53) 



The equivalence is now established by the use of the permutation property given in Eq. (fTTj) . 
In a similar fashion, we find that 

where we must use, see Sect. IIIIBl 



(54) 



dT 



dTo 



0/ V 



dP 



dTo 



dTo)y\dT Jy\dT ,y 



dTo 



dS 



dTo 



dToJy\dT Jy\dV 



dTo 



(55) 



obtained in a similar fashion as Eq. f l53p . 

It is important at this point to relate Cp with Cp^^ and Cy with Cy,^- For this, it is 
convenient to consider the differential dS by treating 5* as a function of Tq, Pq and ^. We 
find that 



Cp = Cp. + T{^] 



d^ 



Cv = Cv, + T[^^ 



Po,V 



dTo 



(56) 



C. Compressibilities Kt and Kg 



The two important isothermal compressibilities are 



Kt = 



dV 



V \dPo 



1 /dV 



V \dPo 



which we need to relate to the corresponding adiabatic compressibility 



dV 



V \dPo 



0/ s 



1 /dV 
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However, we first consider the relation between the compressibility and the expansion coef- 
ficient. We find that 

Kt d{V,To)/d{To,Po) d{V,To) f dTo 



ap d{V,Po)/d{To,Po) d{V,Po) \dPo/y 
The same discussion can be applied to K^^^ and ap^c with a similar result 

Kt,^ _(dT, 



The relation between and -R't,^ and between Ks and Ks^^^ are obtained by treating V 
as a function of Tq, Pq and ^ and of S, Pq and C,, respectively. Using 

dV\ _ fdV\ fdV\ 



we find that 



which is similar to similar relations for the heat capacity in Eq. ( I56|) . We similarly find that 

Let us consider {dV/dPo)rp^ ^: 

dV \ diV, To, d{V, S, d{Po, S, diV, To, fdV\ Cp,^ 



dPoJj.^^^ d{Po,To,0 d{Po,S,Od{Po,To,0 d{V,S,0 \dPo J s,^ Cy,^' 

Similarly, we find that, see Sect. IIIIBl 

dV \ _ fdV\ CV 
dPo)^^~\dPo),C^ 

Thus, we have the standard identity for both kinds of compressibility: 

C p^^ ^T,^ C p Kt 

Let us again consider Ks,^- Rewriting 

(dV\ ^ diV,S,0/diPo,To,0 ^ (dV\ _ idS/dPo)T,,^idV/dTo)p,,^ 
[dPoJs,^ d{Po,S,0/diPo,To,0 ^^0^^,^ idS/dTo)p,,^ 
we find that 



(61) 
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VIII. PRIGOGINE-DEFAY RATIO 



Let us consider Figs. [2] and [3] again that describe various kinds of glass transitions: the 
apparent transitions at Tqq (point D) and Tq^^ (point C) and the conventional transitions 
at Tqq (point D) and Tog (point B). From the discussion in Sect. [TTl we know that the 
Gibbs free energies have a discontinuity between the two states involved at the apparent 
transitions. Even the volumes and the entropies exhibit discontinuities at these transitions. 
On the other hand, the Gibbs free energies, volumes and entropies have no discontinuities 
at the conventional transitions at Tqq and Tog due to the continuity of the state. Let us 
introduce the difference 

Ag = gi - qu (62) 

for any quantity g at a given Tq, Pq in the two possible states I and II. For the apparent 
glass transition at Tqq, qi,qii are the values of q in GL and L, respectively, at Tqq; for the 
apparent glass transition at Tq^\ qj, qu are the values of q in gL and L, respectively at Tq^\ 
For the conventional glass transition at Tqq, qi, qu are the values of q in the glass GL and 
gL, respectively, at Tqq; for the precursory glass transition at Tgg, qi, qu are the values of q 
in gL and L, respectively at Tog. These states are summarized in the Table below. 

Table: Various States 
Apparent Tqq Apparent Tq^-* Conventional Tog Conventional Tqq 

I GL gL gL GL 

II L L L gL 

In terms of the discontinuities ACp, AKt and Aap, the Prigogine-Defay ratio Qj is 



traditionally defined as |Jll3|, l23|-l25( 



ACp AKt 



VTQiAapr 

where it is assumed that the volume is the same in both states at Tq,Pq, as is evident 
from earlier work. As we will see below, the volume is normally not continuous at the 
apparent glass transitions, used in most experimental and theoretical analyses of the glass 
transition. To allow for this possibility, we will consider the following equivalent definition 
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of the Prigogine-Defay ration in this work: 

^ ACpAKt 

~ To{AVap){Aap)' ^ ^ 

where we have absorbed V in one of the Aap-factors. It is clear that 11 is not different from 
-Qtrad ^i^QYL the volume is the same as happens for conventional transitions. 

As the experimentalists have no control over the internal variables, and can only ma- 
nipulate the observables X by controlling the fields yo of the medium, we will discuss the 
evaluation of the Prigogine-Defay ratio in the subspace of yo of the complete thermodynamic 
space of yo, ao. We will consider the simplest possible case in which the subspace reduces to 
the Tq-Po plane. Therefore, we will restrict ourselves to this plane in the following, knowing 
very well that the GL and gL are also determined by the set ^ of internal variables; see Sect. 
IIIIBi We will consider the general case of several internal variables C,k,k = 1, 2, ■ ■ ■ , n. 



A. Conventional Transitions at Tog and Tog 

We will first consider the Prigogine-Defay ratio Ilg at the conventional transitions at 
points B and D (see Figs. [2] and [3]). The continuity of the state across B and D means 
that E, V and S remain continuous across the conventional transitions at B and D. This is 
consistent with the continuity of the Gibbs free energy. Let us first consider the transition 
at B, where the relaxation time r of the system becomes equal to the observation time- 
scale Tobs, so that both states gL and L remain in equilibrium with the medium. Thus, 
T = To, P = Pq, and A = Ao = for both states at B. Therefore, there is no need 
to consider the internal variables in the Gibbs free energy, as they are not independent 
variables. Moreover, V = (9G'/9Po)to and S = —{dG /dTQ)p^y Thus, the Gibbs free energy 
and its derivatives with respect to Tq, Pq are continuous at B; the second derivatives need not 
be. It is clear that B represents a point that resembles a continuous transition in equilibrium; 
it turns into a glass transition curve Tog(Po) of continuous transitions in the Tq -Pq plane. 

For the transition at D, we have a glass GL on the low-temperature side, and gL at the 
high temperature side; both states are out of equilibrium and have the same temperature T{t) 
and pressure P(t), different from To,Po, respectively at the transition. Similarly, A{t) ^ 
is the same in both states. The important characteristics of the conventional transitions 
are the continuity of V and 5 at B and D. We now follow the consequences of these 
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continuities. 



1. Continuity of Volume 

From the continuity of the volume, we have 

dA, \nV = A, ^-^j ^ dn + A, ^-^j ^ dPo = 0, (64) 

where Agq denotes the difference in Eq. (162|) at the conventional glass transitions, and the 
derivatives are also evaluated at the transition points. This equation can be written in terms 
of the compressibilities and the expansion coefficients in the two states at the glass transition 
temperature Tog or Tqq: 



dTo 



dP^ 



^■"^ (65) 



tr Agap 

the isothermal compressibility Kt and the isobaric expansion coefficient ap are given in 
Eqs. ( H5ll and (H71) . respectively, and can be expressed in terms of the derivatives of the 
internal variable ^, such as given in Eqs. ( 159|) and ( 160|) for a single internal variable ^. We 
make no assumption about these ^-derivatives, such as their vanishing or any assumption 
about freezing of ^ at its value at B; indeed, we expect ^ to change continuously over BC. 
Of course, we must remember that ^ is an independent thermodynamic variable in gL and 
GL states only, and not in the L state . The slope equation (l65l) determines the variation 
of Tog or Tog with the medium pressure Pq along the glass transition curve Tog,G(-Po) in 
the To -Po plane, regardless of ^. The form of the above equation does not depend on the 
number of internal variables, provided we use the proper definitions of Kj^ and ap as given 
in Eqs. ( I45l) and f j47l) . respectively. Its form follows from the continuity of the volume at 
the conventional glass transition. 



2. Continuity of Entropy 



From the continuity of the entropy at Tog, we similarly have 



dA,S = Ag ( ^ 



dTo + A„ 



dS 
dPn 



dPn = 0, 



Po V-^O/Tog 

from which we obtain at the precursory glass transition at B 



~dK 



ToAg(V^ap) VgToAgCtp 



tr 



A,C 



gOp 



A.C 



gop 



(66) 



(67) 
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where we have used the equihbrium Maxwell relation {dS/dPo)To = —{,dV/dTQ)p^ = Vap; 
see Eq. ([7]) or Eq. fj4T]) applied to this case. Here Vg is the common volume of gL and L at 
B and has been taken out of Ag{Vap). Again, this relation for the slope is quite general, 
independent of the number of internal variables in gL state at lower temperatures Tq < Tog. 
Accordingly, 

as expected for equilibrium states. It is a consequence of the glass transition being a con- 
tinuous transition between equilibrium states at B. As we will see below, it is not merely a 
consequence of the continuity of volume and entropy simultaneously. 

Let us now consider the glass transition at Tqq. It follows from Eq. fHS]) that 

, rdS\ AgCp 



dToyp, 



In conjunction with Eq. (14T|) . we find that 



dT, 



OG 



VoTAgap {dP/dPo] 



V 



AgCp idT/dTo)s' 

where Vq is the common volume of gL and GL at D and has been taken out of Ag{V ap). 
We finally obtain 

_ AgCpAgKr _ T {dP/dP,)y ^ 
^ VoToiAgapf To {dT/dT,), ^ ^ ' 

for the conventional glass transition at D. The deviation of IIq from unity is independent of 
the number of internal variables. It will be different from unity even if we have no internal 
variables. 



B. Apparent Glass Transitions at Tq^ and Tqg 

Unfortunately, it is not a common practice to determine the Prigogine-Defay ratio at 
the conventional transitions at temperatures Tog(Po) or ^og(-Po)) which resemble continuous 
transition in that the volume and entropy are continuous, along with the Gibbs free energy. 
In experiments, one determines the ratio at apparent glass transitions either at D or at 
Tq^\Pq) in the glass transition region BD; see Figs. |2] and [3l In these transitions, there 
are discontinuities in the G, V and S. The extrapolated point C (see Fig. [2]) identifies 
the apparent glass transition temperature Tq^\Pq), which cannot be treated as a transition 
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temperature because the Gibbs free energy in the two states (gL and L) are not equal, as 
is clearly seen in Fig. |3l It is clear from Fig. |2]that the volume is also different in gL and 
L at the apparent glass transition Tq^\Pq). The discontinuity of the volume should not 
be confused with the continuity of the extrapolated volumes used to determine the location 
of the phenomenological glass transition Tq^\Pq). The extrapolated glass volume does not 
represent the physical volume of the glass at Tq^\Pq) given by the point on the curve BD 
in Fig. [2j The discontinuity is between the physical volumes of gL and L at Tq^\Pq). 
We already know that both the entropy and the enthalpy of the glass continue to decrease 
during vitrification as the system relaxes Indeed, the volume of the glass or gL also 
relaxes towards that of the supercooled liquid L. This will also be true at Tq^\Pq) so that 
the volume and the entropy of gL are higher than their values in the supercooled liquid 
at Tq^^(Po) in a vitrification process. The same sort of discontinuities also occur at D. In 
the following, we will take into account these discontinuities in the volume and entropy 
in determining the Prigogine-Defay ratio at the apparent glass transitions at Tog(-Po) and 
TQg\Po). The discontinuity of volume Ag^^V" (7^ 0) causes a modification of Eq. (IMll at 
these transitions: 

dTo 



(A) , a(A)i^_ a(A).. _ 

1 + 5f Vp) (70) 



7a)~ - 'TiA) 



tr Ag 'ap Ag 'ap 



dPo 
in terms of 

6 In Vj,^^ = rfA(^) In V/dPo \ (71) 

at Tq^"* or Tog, as the case may be; the three Ag'^'''s are the difference A in Eq. ( I62l) evaluated 
at Tq^-* or Tqg, and the new quantity 6g^^Vp has an obvious definition: 

at the appropriate temperature. This contribution would vanish under the approximation 
A'^^ \nV ^ 0, or 6\nVp^^ ~ 0. The slope equation flTOjl must always be satisfied at the 
apparent glass transition temperature. The quantity 6 In V^"^^ in it represents the variation 
of the discontinuity 

Af ) In y = In Vi{To, Pq) - In Vu{To, Pq) 

with pressure along the apparent glass transition curve tI^^\Pq) or Tog(Po), and can also be 
found experimentally. Indeed, we can treat A^^^ In V as a function of Pog = Po(Tq^^) along 
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the transition curves. Then the contribution from the volume discontinuity is given by 



6lnV, 



(A) 



1 dVM 



Vi dPo 



1 dVniPo 



tr 



Vu dPo 



(73) 



tr 



We can use Eqs. (!59l) and (160|) to express the slope in terms of Agi^^.^ and AgOp^^: 



dTo 



dPo 



6 In ri^) + A't'KT,^ - V^,i di/dPo\,, /Vi 



.(A) 



tr 



(74) 



where V^,g represents the derivative {dVi/dS,)rp^^ p^^, and Vi is the GL volume at Tqg or the 
gL volume at T^g^ . The ^-contribution from the L state is absent due to the vanishing of 
the affinity Ao(= 0) in the L. 

Let us now consider the differential of the entropy difference at the apparent glass tran- 
sition in the Tq -Pq plane: 

dS 

Po 



c^A(A)5 = A(A) (^)^^ro + Af) 



dPo 



dPn 



To 



from which we find that 



dTn 



dPn 



54^) - A^-^ (dS/dPo)^^ 



tr 



Af ) {dS/dTo)p^ 



with 



(75) 



(76) 



it represents the rate of variation of the entropy discontinuity 



Af)5 = 5i(To,Po)-5n(To,Po 



along the apparent glass transition curves. Following the steps in deriving Eq. 
that the contribution from the entropy discontinuity is given by 



6S 



(A) 



dPn 



dSujPo 

dPn 



3]), we find 



(77) 



The derivative {dSi/dPo)j,^ in the second term in the numerator in Eq. ( I75|l can be 
manipulated as in Eq. ( |4T1) : 

diS,To) diS,To) diPo,V) fdV\ diS,To) 

diPo,n) 9(r,Po)9(Po,To) [dToJp^^ diV,Po)' 
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in which the last Jacobian reduces to unity under equihbrium by the use of the Maxwell 
relation d{V,P = Pq) = d{S,T = Tq). We, therefore, write 

diS,To) _ {dP/dPo)v _ I 
d{V,Po) {dT/dTo)s 

for the glassy state; this equation also defines the modification 6Syp given by 

I {dP/dPo)v 

- wmh ~ ^^^^ 

for the glassy state, where T, P are the internal temperature, pressure of the glassy state. 
It vanishes under the approximation T = Tq and P = Pq. We now have 

dSi 



(^)p ^ ^^^^^ ^ "^'^^ ^ ^S\,s)a'p. 



For the supercooled liquid, which represents an equilibrium state, we evidently have (^n = 
so that 

Af ) (95/aPo)To = -^i^\Vap) - Vja'p6S'ys- 

We now turn to the denominator in Eq. ( |75l) . For the supercooled liquid state, whose 
temperature is Tq, we have 



we must use Tq^^ or Tqq for Tq to evaluate this slope at the appropriate apparent glass 
transition. For the glass, whose internal temperature is T, we have 



where we have introduced a correction parameter 

6T'^T^^},/T-1, (79) 

with Tq^q denoting Tq^^ or Tqq as the case may be. Again, this modification term vanishes 
under the approximation T = Tq. We thus find that 

ToAf ) {dS/dTQ)p^ = Af )Cp + C'pdT'. 
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Equating the two different versions of the slope in Eqs. (170|) and (1750 . we have 



[1 + 6'U,a), 



where we have introduced a new quantity ^'IlgA, whose definition is obvious from the equahty. 
We finally find that the Prigogine-Defay ratio is given by 

Tg^g^^jA^ 'apA^ '(Vap) 1 + 5^ 'Vp 

at the apparent glass transition. Its complete form is given by 

1 + (Via'pSS'ys + ^4^^) /Ai''\Vap) 



It should be obvious that the Prigogine-Defay ratio is itself a function of time as it depends 
on time- dependent quantities such as Ag^''5', etc. 



1. Approximation A 

Let us assume that the discontinuities in the volume and entropy are negligible or that 
the contributions 6 In Vp^'' and 6S^p^^ are negligible. In that case, the Prigogine-Defay ratio 
reduces to 

l + V,aUSl.s/Ai^\Vap) 
1 + Cj^ST^Ai^^Cp 

and will have a value different than 1. Thus, the continuity of volume and entropy alone 
is not sufficient to yield IlgA = 1, as noted above. If we further approximate T ^ Tg 
and P Pq, then SSyg and 6T^ 0, and we obtain IlgA — 1- This is expected as 
the approximations change the apparent glass transition into a continuous transition. If, 
however, we only assume P ~ Pq? but allow T to be different from Tq, then 

1 

^^""^ - {dT/dT,)s, 

and we still have IlgA 7^ 1- 
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2. Approximation B 



We make no assumption about 5 In V^^^ and 5S^p \ but approximate T ~ Tq and P c:^ Pq. 

obtain 



In this case, 5Syg ~ and 5T^ ~ 0, and we obtain 



1 + 5f Vp ■ 

If, however, the approximation T Tq is not valid, we have 

iigA — 

(1 + ^ 

In both cases, IlgA 7^ 1. 

C. Comparison with Other Attempts for 11 



~ (1 + Cl5TG//\[^^Cp){l + 5f Vp) 



As far as we know, almost all previous attempts |9l-ll3l. |23-25| in the evaluation of 11 are 
based on treating the glass transition as a direct transition from L to GL; the structure is 
supposed to be almost frozen in the latter. As we see from Figs. |2]and|3l this can only occur 
at C between L and the extrapolated branch DC. At C, there will be a discontinuity between 
the values of the internal variable ^; it will take the equilibrium value ^'q in L, but will take a 
non-equilibrium value ^q^^ 7^ S,'q obtained along DC. Similarly, A = = in L at C, while 
A = Ac 7^ in the extrapolated GL at C. As C is obtained by matching the volumes, the 
volume remains continuous, but there is no reason to believe that the entropy will remain 
continuous in this transition. The Gibbs free energy obviously remains discontinuous in this 
transition. 

However, we have been careful in not treating this transition as an apparent transition 
above for the simple reason that there is no guarantee that the branch DC can be described 
by vitrification thermodynamics at the constant cooling rate r. To see it most easily, we 
observe that as the cooling rate is gradually taken to be slower and slower, the transition 
point B gradually moves towards C along BF. However, the analog of BD will most certainly 
not be identical to DC for the simple reason that the state of L will continuously change 
to gL so that the values of ^ and A in gL at C will be identical to their values ^ = 
and A = in L at C. Moreover, there is no guarantee that the extrapolated branch DC can 
even satisfy thermodynamics with known controllable parameters Tq, Pq and r. To treat this 
"transition" as a glass transition requires some approximation, which we have avoided. 
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IX. CONCLUSIONS 



We have followed the consequences of internal equilibrium to derive generalizations of 
equilibrium thermodynamic relations such as Maxwell's relations, Clausius-Clapeyron re- 
lation, relations between response functions (heat capacities, compressibilities, etc.) to 
non-equilibrium systems. Non-equilibrium states are described not only by internal fields 
(temperature, pressure, etc.) that are different from the medium, but also described by in- 
ternal variables which cannot be controlled from outside by the observer. The observer can 
only control the observables. Thus, in this work, we have also discussed how the thermody- 
namics should be described in the subspace of the observables only. As glasses are a prime 
example of non-equilibrium states, we have reviewed the notion of the glass transition. The 
frozen structure known as the glass (GL) does not emerge directly out of the equilibrium 
supercooled liquid (L). There is an intermediate non-equilibrium state (gL) that is not yet 
frozen when it emerges continuously out of the equilibrium liquid L. At a lower temperature, 
this state continuously turns into GL. Because of this, we find that there is no one unique 
non-equilibrium transition. We introduce four of the most conceptually useful transitions. 
At two of them, which we term conventional glass transitions, the Gibbs free energies and 
the states are continuous. Thus, they are the non-equilibrium analog of the conventional 
continuous or second order transition between equilibrium states. At the other two glass 
transitions, which we term apparent glass transition, not only the states but also the Gibbs 
free energies are discontinuous. Because of this, these transitions are examples of a zeroth 
order transition where the free energy is discontinuous. But there is no transition in the 
system itself at the apparent glass transition as discussed in Sect. [ITl 

We briefly review the use of Jacobians which are found extremely useful in obtaining 
the generalization of the Maxwell relations. There are many other Maxwell relations than 
reported here; they can be easily constructed. We then discuss various response functions 
and obtain relationship between them in non-equilibriums states. Surprisingly, many of 
these relations look similar in form to those found in equilibrium thermodynamics. 

We finally evaluate the Prigogine-Defay ratio at the four possible glass transitions. We 
find that the ratio is normally different than 1, except at the conventional glass transition at 
the highest temperature, where it is always equal to 1, regardless of the number of internal 
variables. We also find that the continuity of volume and entropy is not a guarantee for 
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n = 1. We compare our analysis of H with those carried out by other workers. 
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Appendix A: Relation between 2- and 3-Jacobians 



Let us consider a function F{x,y,z), where x,y,z may stand for To,Pq,^, respectively. 
Then 

dF — F^y^dx + Fy^zj-dy + F^^^ydz, 
where we have used the compact notation 

'OF' 



F = 

^ x,yz — 



dx 



yz 



etc. Now, 



Similarly, 



dF\ 



( dz 



\dxj' 



j ~ ^x,yz + Fz^xy I 

OK 



dy 



K +K [-] 

^^y,zx ^ ^^z,xy \ Qy J ' 



We express Fx^yz as a 3-Jacobian and manipulate it as follows: 



d{F,y,z) ^ d{F,y,z) d{K,x,z) ^ d{F,y,z) 
d{x,y,z) d{K,x, z) d{x,y, z) d{K,x,z) 



dK 
dy 



z,xy 



dz 
dy 



dK 



Jx 



d{F,y,z) 
d{K, X, z) 



K 



z,xy 



dz \ d{F, y, z) 



dKj^d{K,x,z) 



Using this, we find that 
'dK 



dF\ 



dy 



d{F,y,z) 
d{K, X, z) 



K 



z,xy 



dz \ d{F, y,z) d{z, y) 

+ i'z 



dKj^d(K,x,z) 



z,xy 



d(K, x) 



Let us call the quantity in the square brackets D, which can be rewritten as 

d(F. y) , 



where 



d{K, x) 



j^, ^ d{F\y^d{K^ _ d{z,x) d{F,y,z) d{z,y) 
d{K,x,z) d{F, y) ^'^^ a(F, y) d{K, x, z) ^ ^'^^ d{F, y) ' 
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Using Eq. (fl^ . it can now be shown in a straight-forward manner that 



which proves that 

(^] =-(—) ^^^'^^ 
\dxjy \dyJ^d{K,x)' 

the desired result. 

While we considered F and K as a function of 3 variables, we can generalize the result 
to any number of variables. We will not pause here to do that. 
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